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It is widely accepted that the broad absoption hne region (BALR) exists in 
^ ■ most (if not all) quasars with a small covering factor. Recent works showed that 

^ . the BALR is optically thick to soft and even medium energy X-rays, with a typical 

O ' hydrogen column density of a few 10^^ to > 10^*^ cm~^. The electronic scattering 

^ . in the thick absorber might contribute significantly to the observed continuum 

. polarization for both BAL QSOs and non-BAL QSOs. In this paper, we present a 

^ . detailed study of the electronic scattering in the BALR by assuming an equatorial 

p . and axisymmetric outflow model. Monte-Carlo simulations are performed to 

Q . correct the effect of and radiation transfer and attenuation. Assuming an average 

+i , covering factor of 0.2 of the BALR, which is consistent with observations, we find 

. the electronic scattering in the BALR with a column density of ~ 4 x 10^^ cm~^ 

> . can successfully produce the observed average cotinuum polarization for both 

^ , BAL QSOs and non-BAL QSOs. The observed distribution of the continuum 

polarization of radio quiet quasars (for both BAL QSOs and non-BAL QSOs) is 
helpful to study the dispersal distribution of the BALR. We find that, to match 
the observations, the maximum continuum polarization produced by the BALR 
(while viewed edge-on) peaks at P = 0.34%, which is much smaller than the 
average continuum polarization of BAL QSOs (P = 0.93%). The discrepancy can 
be explained by a selection bias, that the BAL with larger covering factor, and 
thus producing larger continuum polarization, is more likely to be detected. A 
larger sample of radio quiet quasars with accurate measurement of the continuum 
polarization will help give better constraints to the distribution of the BALR 
properties. 



Subject headings: polarization-scattering-quasars: absorption lines 



- 2 - 



1. Introduction 

About 10-20 % optically selected QSOs exhibit broad absorption troughs in resonant 
lines up to 0.1c blueward of the corresponding emission lines (Hewett & Foltz 2003, Reichard 
et al. 2003). Usually the Broad Absorption Lines (BAL) are detected only in high ionization 
ones, such as CIV, NV, SilV and OVI, named high ionized BAL (HiBAL), but 10% of 
BAL QSOs show also low ionization lines (LiBAL), such as Mgll, AlIII and even FelL The 
blue shift of the absorption lines suggests that they are formed in a partially ionized wind, 
outflowing from the quasar. The observed flux ratios of the emission to absorption line imply 
that the covering factor of the BAL Region (BALR) be < 20% (Hamann, Korista & Morris 
1993). Futhermore, the properties of UV/X-ray continuum and emission lines of BAL QSOs 
and non-BAL QSOs are also found to be similar (Weymann et al. 1991; Reichard et al. 2003; 
Green et al. 2001). These facts lead to a general picture that BALR covers only a small 
fraction of sky and may be present in every quasar (Weymann et al. 1991). We note there 
are evidences suggesting that LiBAL QSOs are in a special evolution phase of quasars rather 
than merely viewed on a special inclination (Voit, Weymann & Korista 1993; Canahzo & 
Stockton 2001; c.f., Willott, Rawlings,& Grimes 2003; Lewis, Chapman, & Kuncic 2003). 

The evidence for non-spherical BALR is also supported by the discovery that BAL 
QSOs are usually more polarized than non-BAL QSOs (Ogle et al. 1999). BAL QSO is the 
only high polarization population among the radio quiet QSOs (e.g.. Stockman, Moore & 
Angel 1984). Stockman et al. found 9 of 30 BAL QSOs show high polarization [P > 1.5%) 
in the optical continuum. The result was conflrmed later by Schmidt & Hines (1999), who 
found the average polarization degree of BAL QSOs is 2.4 times that of the optically selected 
quasars. Consistently, Hutsemekers & Lamy (2001) obtained average polarization degrees of 
~ 0.43%, 0.93% and 1.46% for non-BAL QSOs, HiBAL QSOs and LiBAL QSOs respectively. 
This result is similar to that of Schmidt & Hines (1999): 0.4% and 1.0% for non-BAL and 
BAL QSOs. It is also evident that the polarization distribution of non-BAL QSOs drops 
sharply toward high polarization (see also Berriman et al. 1990). 

BAL QSOs are notorious X-ray weak following the work of Green & Mathur (1996; see 
also Brinkmann et al. 1999). Now we have strong evidences that the weakness in X-rays is 
not intrinsic but due to strong X-ray absorption, with the hydrogen column density of a few 
10^3 to > 10^4 cm-2 (Wang et al. 1999; Gallagher et al. 2000; Mathur et al. 2000). The X-ray 
absorber might be responsible for the recently detected blueshifted X-ray BALs, suggesting 
they are outflowing at even higher velocities and highly ionized (Chartas et al. 2002; 2003). 
The electronic scattering in the thick X-ray absorber may contribute signiflcantly to the 
observed continuum polarization for both BAL QSOs and non-BAL QOSs, if BALR exists 
in every quasar. In this paper, we perform detailed calculations and Monte Carlo simulations 



-3- 



to study the polarization produced by electronic scattering in the BALR. By comparing the 
expected polarizations with the observed ones, we give strong contraints to the BALR model. 

2. Models and the Monte-Carlo method 

There are many dynamic models for the BAL outflow, depending on the flow type 
(hydrodynamic or hydromagnetic flow) and the accelerating mechanisms (radiation, gas 
pressure or magnetic field). In most models, the flow is accelerated through the resonant 
line scattering, which is supported by the line-locking phenomena. A general difficulty of such 
models, however, is to prevent the flow from being fully ionized while its average density drops 
rapidly with increasing radius. One solution proposed by Murry & Chiang (1995, hereafter 
MC95) is that the flow is shielded from the intense soft X-ray radiation by highly ionized 
medium in the inner region with a typical column density of a few 10^^ cm~^ to 10^^ cm~^. 
This shielding gas can also account for the observed heavy X-ray absorption. This model is in 
qualitative agreement with more recent hydrodynamic calculation of the radiative accelerated 
wind from an accretion disk (Proga, Stone & Kallman. 2001). The second solution is the 
two-phase flow, in which a dense, low ionization, cold clouds are embeded in a highly ionized, 
hot and tenuous medium. The cold clouds with a small filling factor, accelerated by the line 
and continuum radiation pressure to high velocities, is responsible for the BAL features. 
An implement to the second scheme is the massive hydromagnetic and radiative driven wind 
model. In the model, the massive high-ionized continuous outflow is driven centrifugally, and 
accelerated radiatively by the central continuum source (Everett 2002; Konigl & Kartje 1994; 
de Kool & Begelman 1995). The total column density of the hot medium is also very large 
(with Nh — 10^^ cm~^ ~ 10^^ cm~^ ). In fact such a medium itself can also be considered 
as the shielding gas. A third scheme is the dusty wind due to the mass loss of stars in the 
nucleus; both the dust absorption and the line scattering contribute to the accelerating force 
of the gas (Voit et al. 1993; Scoville & Norman 1995). 

Following MC95, we assume an equatorial and axisymmetric outflow, with a half open 
angle 9o (see Fig. 1). In this paper, we focus on the continuum polarization produced by the 
electronic scattering in the outflow, and leave the study of the resonant scattering, which 
can produce obvious polarization around the broad absorption trough, in a future paper 
(Wang et al. in prep). The shielding gas is the major source of the electronic scattering 
and X-ray photo-electronic absorption. Since both the scattering and X-ray absorption are 
insensitive to the density profile of the shielding gas, we assume a constant electron density 
in this region. We consider a range of column densities (10^^ to > 10^^ cm~^) for this region, 
consistent with that obtained from the X-ray observations. 
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FoUowing Lee, Blandford & Western (1994), we denote the density of the polarized 
incident photons in direction {9i,(fii) as 

V P21 P22 / 

where pij is the i,j component of the photon-density matrix. The outward density in 
the direction {9o, (po) is 

( Pii P12 \ 

V P21 P22 J ' 

For the electronic scattering we may write the outward photon density after one scat- 
tering as (see Chandrasekar 1950): 

pIi oc p\i [cos^ 9o cos^ 9i cos^ {(pi — (po) + 2 cos 9o cos 9i sin 9o sin 9i cos{(pi — (po) + sin^ 9o sin^ 9i 

+PI2 cos^ 9o sin^((pi - ipo) + ^{Pi2 + P21) [cos^ cos 9i sm{2(pi - 2(po) 
+2 cos 9c, sin 9o sin 9i sin((^j — <^o)] 

P22 oc p\i cos^ 6'j sin^((/7i - (po) + P22 cos^((/7j - (/?„) - ^(^^2 + P21) cos 6lj sin(2(/7j - 2(po) (2) 
P12 Pii [ — cos 9o cos^ 9i sin(2<^j — 2(po) — sin 9o sin cos 9i sin(y?.j — (po)] 

1 ■ 1 1 

-I — P22 cos 9o sin(2(/7j — 2(po) + pli[—- cos cos 9i + - cos cos 9i cos(2(/7j — 2</7o)] 
. 1 1 

+Pi2[~ cos 6'o cos6'i H — cos 9o cos 9i cos(293j — 2990) + sin^^o sin 9i cos{pi — Po)] (3) 
2 2 

P21 oc p*^]^ [ — COS 9o cos^ sin(2(^j — 2po) — sin sin 9i cos sin((^j — po)] 
2 

+ ~P22 cos 6*0 sin ( 293 j — 2po) + Pi2[ — cos 9o cos 6*1 -I — cos^^o cos^i cos{2pi — 2po)] 

2 2 2 

+P2i[^ COS 9o COS 9i + l- COS COS 9i cos(2</7j — 2</7o) + sin 9o sin cos(</7i — po)] (4) 
The four STOKES parameters read. 



1) 



I ^ P11+ P22, <5 = Pii - P22, = P12 + P21, V ^ pi2 - P21 (5) 
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Other interesting quantities, such as the total flux, the PA rotation, the polarization 
degree and the polarized flux, can be calculated from the Stokes parameters. The polarization 
degree follows 



P = J (6) 

Following the steps below, we run Monte-Carlo simulations to calculate the output 
Stokes parameters for given incident radiation and spatial distribution of the scatterer. 



• For each incident photon, we give a random propagation direction. Then, the unit 
density matrix of the photon with a given polarization degree is calculated (Lee 1994). 

• The probability of the photon passing through an absorber with an optical depth 
T is p = exp(— r). By assigning a random number in [0,1] to p, we calculate the 

corresponding r. If r is larger than the real optical depth along the photon propagation 
direction (r'), the photon will escape from the medium and be collected in the output 
basket. 

• If r' > T, the photon will be scattered at r, where the real optical depth along the path 
of the photon is r. The emergent direction and the density matrix of the scattered 
photon is calculated (equations 1,2,3 &4). 

• The frequency of the scattered photon is calculated by taking into consideration of 
Doppler shift, which is a function of the incident direction, emergent direction and the 
velocity vector of the scattering particle. The velocity vector of the scattering particle 
is determined by the bulk velocity and the thermal velocities at r (Lee & Blandford 
1997). 



We repeat step 2, 3, 4 until the photon either escapes or is absorbed by the accretion 
disk, which is assumed to be optically thick with no reflection. 



3. Results and Discussion 

In the case of optical thin limit (r << 1), the polarization degree viewed at an inclination 
angle i can be written analytically as (Brown & McLean 1978, hereafter BM78) 



Pe = 2f(l - 37) sm^i 



(7) 
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where 



111 r{fJ')df^ 



7- / . (8) 



and 



(9) 



where ^ = cos 9 and t(//) is the Thomson depth at different polar angle 9. 

If BALR covers the inclination angle from 90° - 9o to 90° + 9o (see Fig. 1) and the 
accretion disk is optically thick with no reflection, the average polarization due to the electron 
scattering for BAL and non-BAL QSO can be written as 



<pnon-BAL^ 2/3 - //q + //q/S 

where /xq = sin 6*0- 

Assuming a constant Thomson depth (tq) of the the shielding gas at different directions, 
we obtain the average polarization for BAL QSOs 

^B = (l-^/X^)^To//o(l-/^o) (11) 

The constant column density of the shielding gas at different direction is over-simplified. 
However, a distribution of r will not significantly change the estimated mean optical depth 
{J^"^^T{fi)dfi) if sin^o is small. 

Using the observed values of R = 2.2 and Pb = 0.93% (see §1), we obtain /xo=0.5 and 
ro=0.144. /io = 0.5 is much larger than 0.2 derived from the observed fraction of BAL 
QSOs. We point out the discrepancy might be due to selection bias, if iiq has a dispersal 
distribution instead of a 5 function for quasars. We can see that we have higher chance to 
detect BAL QSOs with higher //q, and reversely, higher chance for non BAL QSOs with lower 
Hq. Thus the average no of the observed BAL QSOs tends to be higher than that derived 
from the fraction of BAL QSOs (also see Morris 1988), and that of the observed non-BAL 
QSOs tends to be lower. For example, an average /Iq — 0.25 for the observed BAL QSOs 
and /Iq — 0.188 for the observed non-BAL QSOs can easily macth the observed value R = 
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2.2. Using Ho — 0.2, to produce the observed average polarization of BAL QSOs, tq = 0.26 
or Nf> = 3.9 X 10^^ cm~^ is required. If assuming an average fiQ = 0.25 for the observed 
BAL QSOs, we obtained a shghtly lower N,. = 3.25 x 10^^ cm~^. The predicted average 
column density is consistent with those derived from the X-ray observations (Wang et al. 
1999; Gallagher et al. 2000; Mathur et al. 2000). 

Following the procedures described in §2, wc perform Monte-Carlo simulations to correct 
the effect of radiative transfer. Wc test the Monte-Carlo code by simulating single scattering 
process, and found the results from the simulations are in good agreement with the analytic 
ones (Fig. 2). In the simulations below, wc consider the half covering angle from 6° to 
40° and the column density of electron from 10^^ cm~^ to 7 x 10^^ cm~^. In Fig. 3 we 
plot the resulted average continuum polarization degree (Pb) for BAL QSO as a function of 
the column density. We can see that, to produce the observed average polarization degrees 
of HiBAL (0.93%), the minimum column density of 2 xlO^^ cm~^ is required for the range 
of 6o considered. This number increases to 3.5 xlO^^ cm~^ for LiBAL. Using the average 
9o = 12° (corresponding to /Xo=0.2), we derive a column density of 4 xlO^^ cm~^ for HiBAL, 
slightly higher than the analytical result. In Fig. 4, we plot the output polarization degree 
as a function of the viewing angle from Monte Carlo simulations, assuming A^^g = 4 xlO^^ 
cm~^ (re = 0.266). For comparison, the analytic results, which is valid when Tg << 1, are 
also presented. The impact of the radiative transfer to the polarization degree is apparent 
in the figure. Basically, because of the attenuation of the direct continuum, the radiative 
thransfer causes larger polariztion for BAL QSOs than that from the analytic calculation. 
For non-BAL QSOs, the situation is more complex. For a model with small Oq — 12°, the 
optical depth at small inclination angle is smaller, consequently the scattered photons escape 
more likely along the polar direction, which results in slightly larger (smaller) polariztion at 
smaller (larger) inclinations. However, for a model with larger 6q = 29°, the attenuation of 
the scattered light becomes important, which reduces the polarization for non-BAL QSOs. 

We point out that, in additional to the observed average polarization degrees for BAL 
and non-BAL QSOs, the distribution of the observed polarization of a radio quiet QSO 
sample can provide further constraints to the model of the scatterer. We consider the 
optically thin case first. For a given axisymmetric electron-scatterer, the polarization degrees 
at different view angles i follow eq. 7. The normalized distribution of the polarization degree 
can thus be written as 



dN 1 1 ^ 

—r- = , / , = , for p <po (12) 

dp dp/dji 2poVl -p/po 

where po is the polarization degree viewed at the direction perpendicular to the sym- 
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metric axis. Note po depends on the optical depth and the geometry of the scatterer. We 
notice that the distribution is a monotonic increasing function of p for a given scatterer. 
Monte-Carlo simulations, which taking account of radiative transfer for larger optical depth, 
give similar result (see. Fig. 5). However, an opposite trend was found in the observed 
distributions given by Stockman et al (1984) and Berriman et al. (1990), who found the 
observed distribution of polarization degree of all QSOs (including BAL QSOs) decreases 
with increasing p above 0.2%. We argue that the discrepancy might also be explained by a 
dispersal distribution of the properties of the shielding gas. Here considering a distribution 
of Po as fipo), we can write formally the distribution of p as 

dN^ fipo)dpo 
dp Jp 2poy/l -p/po 

If the observed distribution of the polarization is known, we can solve the above equation 
inversely, e.g., through Richardson-Lucy approach, to obtain the f{po), which in turn can be 
used to constrain the geometry of the scatterer. Using the polarization distribution of QSOs 
presented by Stockman et al (1984, see Fig. 1 in their paper), which can be described 
approximately with 

d(V _ r 16.5 (p < 0.2%) 

'dp~\ 22.23 exp(-p/0.0067) (p > 0.2%) 

we obtain f{po) in Fig. 6, which rises steeply toward lower po with a peak at po — 
0.34%. This suggests that the shielding gas in a large fraction of QSOs either have smaller 
column density, or smaller covering factor. Note that P — 0.34% is much smaller than 
the average continuum polarization of BAL QSOs which is 0.93%. This suggests that most 
of the quasars have BALR with covering factor much smaller than 0.2, thus these BALR 
make less contribution to the observed sample of BAL QSOs because of the smaller chance 
to be detected. A larger sample of radio quiet quasars with accurate measurement of the 
continuum polarization will help give better constraints to the distribution of the BALR 
properties. Note the observed distribution at the small polarization end {p < 0.2%) was 
not well constrained because of the limited accuracy of the measurements, thus the resulted 
shape of f{po) has large uncertainty. 

Note that the derived f{po) is also valid for axisymmctric scattering models, including 
the polar scattering models (e.g.. Ogle 1997; Lamy & Hutscmckcrs 2004) and possible bi- 
conical scattering model (Elvis 2000). The reason is that the polarization produced by 
axisymmetric outflow is cx sin^i, but independent to the outflow model (eq.7). 



(13) 



-9- 



In passing, we point out that the average level, as well as the angle-dependence, of the 
polarization degree may be affected by the anisotropic emission of the continuum source, 
which is neglected in the above calculations. Certain levels of anisotropy are expected in 
optically thick accretion disk models, which usually predict stronger emission in polar direc- 
tions than on the equatorial plane. In such models, the average polarization degree would 
be lowered in our equatorial scatterer model, and the polarization degree increases with in- 
clination faster than in the case of isotropic continuum emission. However, we believe that 
the anisotropic emission would not severely affect our results for two reasons. First, BAL 
QSOs are at least as luminous as non-BAL QSOs ^. Second, we do not find anti-correlation 
between the polarization degree and the continuum luminosity (Moore, R. L. & Stockman, 
H. S. 1984), which is expected in a model with anisotropic continuum. 



4. Conclusion 

We show that the X-ray absorber in the BAL QSOs is capable to reproduce the observed 
continuum polarization for both BAL QSOs and non-BAL QSOs. However, the covering 
factor of the BALR in quasars is required to have a dispersal distribution, instead of a 5 
function. We also find that to macth the observed distribution of the continuum polarization 
of radio quiet quasars, the BALR in most QSOs produces much smaller maximum continuum 
polarization pO (while viewed edge-on) with a peak at pO = 0.34%, which is much smaller 
than the average continuum polarization of BAL QSOs, which is 0.93%. Consequently, the 
BAL QSOs with small pO are likely to have covering factor of BALR much smaller than 0.2, 
thus make less contribution to the observed sample of BAL QSOs because of the smaller 
chance to be detected. 
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Fig. 1. — The geometry of the outflow model used in the paper. 
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Fig. 2. — The output polarization of the single-scattered hght viewed at different inchnation 
angle i in the case of optically thin limit (r << 1). The geometry of the scatterer is presented 
in Fig. 1. Here fj, = cosi, = 12°, and iVe = 4 x 10^^ cm~^. The analytic result is presented 
as open circles (eq. 7), and the simulation resuls as black squares. 
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Fig. 3. — The average continuum polarization degree for BAL QSO as a function of the 
column density for different ^^o^ based on Monte-Carlo simulations. The two horizontal lines 
mark the polarization degree 0.93% (the average value of HiBAL QSOs) and 1.46% (the 
average of LiBAL QSOs). The data 
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Fig. 4. — The effect of the radiation transfer on the output polarization for iVg = 4 x 10^^ 
cm~^ (re = 0.266). The analytic result takes into account of only single scattering while both 
multiple scattering and the attenuation of the direct light is considered in the Monte-Carlo 
simulation. 
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Fig. 5. — The distributions of the polarization degree for models with 6q — 12° and 60 — 29°. 
The left two panels are the analytic results and the two right panels are the simulation 
results. A column desity Ne — 4: x 10^^ cm~^ is assumed. 
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Fig. 6. — The reconstructed po distribution assuming an axisymmetric electron-scattering 
model, where po is the polarization degree detected on equatorial plane. See text for details. 



